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Abstract 

This article shows a strong averaging principle for diffusions driven by discontinuous heavy¬ 
tailed Levy noise, which are invariant on the compact horizontal leaves of a foliated manifold 
subject to small transversal random perturbations. We extend a result for such diffusions with 
exponential moments and bounded, deterministic perturbations to diffusions with polynomial 
moments of order p ^ 2, perturbed by deterministic and stochastic integrals with unbounded 
coefhcients and polynomial moments. The main argument relies on a result of the dynamical 
system for each individual jump increments of the corresponding canonical Marcus equation. 

The example of Levy rotations on the unit circle subject to perturbations by a planar Levy- 
Ornstein-Uhlenbeck process is carried out in detail. 

Keywords: Markov processes on manifolds; solutions of stochastic differential equations with Levy 
noise, foliated manifolds; strong averaging principle; scale separation; Marcus canonical equation; 
dynamical systems; heavy tail distributions; 
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1 Introduction 

The theory of averaging (deterministic) ordinary differential equations, whose origins date back to 
the works of Laplace and Lagrange, has been applied through its history in many fields of appli¬ 
cations such as celestial mechanics, nonlinear mechanics, oscillation theory and radiophysics. First 
rigorous results start with the foundational contributions of Krylov, Bogoliubov and Mitropolskii 
[n la la [371 m]. For a comprehensive and systematic introduction to the subject we refer to 
the monograph [55] of Saunders, Verhulst and Murdock. The idea of averaging random systems 
given as stochastic ordinary differential equations with respect to Gaussian processes goes back to 
Stratonovich [57] (58] , first rigorous results appear with the seminal works of Khasminski and others 
[3I1132113311S1126112I1156] . More recent developments on systems of stochastic (partial) differential 
equations with continuous Gaussian noises can be found for instance in [H [iTl [35l SZl EH [62l (63] 
and the references therein. 

However, in many contexts the Gaussian paradigm is known to be too limited. First results on 
averaging differential equations with respect to discontinuous and non-Gaussian Poisson noise are 
obtained in [6] and [36|. A first strong averaging principle for scalar Levy diffusions with Lipschitz 
coefficients and bounded jumps is established in m- In [28] the authors show a strong averaging 
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principle for Levy diffusions with exponential moments on foliated manifolds, explained in more 
detail below. This article is dedicated to the averaging of a large class of discontinuous semimartin¬ 
gales supported on the compact leaves of a foliated manifold having polynomial moments. 

Intuitively speaking a foliated manifold is a Riemannian manifold equipped with a family of 
equivalence classes of submanifolds depending on a parameter, which defines a “transversal” com¬ 
ponent. Precise details on foliated manifolds are found in classical texts on the topic such as 
[Ml [591 [64]. The notion of a foliated Brownian motion is introduced in the seminal article [24| 
by Garnett. In |33| Li shows an averaging principle for Hamiltonian systems, which inspired first 
results on averaging on foliated spaces in [25], where the authors show an averaging principle for a 
foliated Brownian diffusion. More precisely, the authors study a Brownian diffusion in Stratonovich 
sense on a foliated manifold, which respects the foliated structure of the manifold in the sense that 
the diffusion does not leave the compact leaf of its initial condition almost surely. By assumption, 
the diffusion enjoys a unique invariant measure supported on the leaf of its initial condition as the 
limit of its time average in sense for some p ^ 2. For small e the law of the perturbed diffusion 
converges on the accelerated time scale t/e to the invariant measure, such that the transversal 
perturbation converges to the vector field averaged against the invariant measure, which leads to 
an ordinary differential equation in transversal direction to the leaves of the foliation. Their main 
result is an averaging principle in the sense with logarithmic rates of convergence. 

This result is extended in [28] to a class of foliated Levy diffusions with compact leaves and 
exponential moments of the underlying jump Levy process, formulated in terms of canonical Marcus 
equations, as for instance in [H]. However by the Levy-Chinchine formula (see for instance m) 
it is obvious that this is a rather narrow subclass of possible Levy drivers. The current article 
generalizes this result to the case of Levy jump diffusions, with moments of order p, p ^ 2, for 
which the averaging converges in sense. The lower bound p ^ 2 seems natural for equations on 
manifolds since the Marcus canonical integral can be rephrased as an integral against the quadratic 
variation of the underlying process, see [4T], Lemma 2.1. 

The difficulty of an immediate extension of the results to Levy diffusions with only p-th mo¬ 
ments lies in their formulation as a canonical Marcus equation [msasH], where each single jump 
increment of the process is given as the solution of an ode, with a vector field tangential to the 
leaf. Assume that AZ is a single jump increment in an appropriate noise space of the driving Levy 
process Z. The jump increment on the manifold then has to follow the local coordinates and is 
then given as the increment — y, where = Y(l;y, F AZ) is the time 1 map of 

the solution Y of the ordinary differential equation 

dY 

— = F{Y)AZ, Y{0) = y, (1) 

and T is a Lipschitz vector field such that FAZ takes values in the tangent space of our manifold. 
For details we refer to [Mia HI]. The problem is that in general the Lipschitz continuity of F 
only implies 

^ - y| Vx,y, (2) 

where i is the Lipschitz constant of F. This means each jump increment on the manifold depends 
exponentially on the random size of AZ. See m Lemma 3.1. Taking the expectation of ([2|) the 
hniteness of the right-hand side implies the exponential integrability < oo 

of the Levy measure v of Z for some constant «: > 0 larger than the Lipschitz constant i of F, 
which is equivalent to the existence of exponential moments of Z. This straight-forward argument 
is the main reasoning concerning the moments carried out in the previous article [28]. However, 
since the leaves of the foliation are compact and the main driving diffusion X is invariant on the 
leaf of its initial condition, any jump increment of X is bounded by the diameter of the leaf in 


2 


the surrounding space. With this intuition in mind we may prove in Lemma 13.11 of Section 3 the 
following result on positive invariant ODE dynamical systems V of type dll), which yields a global 
constant O > 0 such that for any x, y 

sup|(DF(F(t;x)AZ)F(y(t;x))AZ- (DF(y(t;y)AZ)F(y(t;y))AZ| ^ C\x - y\ ||AZf. (3) 
t^o 

Taking the expectation of ([3|) a hnite right-hand side is equivalent to /|| 2 ||>i < oo imposing 

only second moments of Z. It turns out eventually to be an easy task to link ([2]) to ([3]) via Taylor 
expansion of We follow these lines of reasoning in a technically more subtle setting in 

Section [3l 

The second extension we undertake is the step from perturbations by a small deterministic 
bounded vector field in [28] to a general class of discontinuous Levy diffusions with moments of 
order 2p, whose multiplicative coefficients may depend on the slow component. The coefficients in 
front of odB and oZ will be only depend on since it is well-known in averaging theory that 

in general diffusion coefficients are difficult to average in a strong sense. 

The article is organized as follows. Subsection 2.1 lays out the general setup. Subsection 2.2 
states the specific hypotheses on the integrability and ergodicity conditions of the stochastic pro¬ 
cesses and the main result of this article given in Theorem 12.21 Subsection 2.3 spells out the 
main example: Levy processes with polynomial moments on the unit circle. In Section 3 we es¬ 
tablish the estimate ([3]) and derive the crucial estimates on the deviation of the perturbed from 
the unperturbed solution, under arbitrary Lipschitz functions including a crucial dynamical system 
argument. Section 4 is dedicated to the control of the averaging error term exploiting the results 
from Section 3 in special cases. Section 5 finishes the proof of the main result synthesizing Section 3 
and 4. The article finishes with an Appendix providing the missing details of the example. 

2 Object of study and main results 

2.1 The set up 

The geometry: Let M be a hnite dimensional connected, smooth Riemannian manifold. It 

is known by the strong version of Whitney theorem for instance in Boothby |I0] that any hnite 
dimensional smooth manifold is embedded in for some m £ N sufficiently large. The manifold 
M is equipped with an n-dimensional foliation 5JI in the following sense. Let 5JI = {Lx)x£M, with 
M = UxGAf equivalence classes of the elements of M satisfying the following 

properties. 

a) Given an xq £ M, there exists a neighborhood [/ C M of the corresponding leaf a 
connected open set R C containing the origin 0 £ and a diffeomorphic coordinate map 
if : U ^ Lxq X V. 

b) The set U of item a) can be taken small enough such that the derivatives of the coordinate 
map (/? are bounded. The second coordinate of a point x £ C/, called the vertical coordinate, 
is denoted with the help of the projection tt : C/ —>• R by ip{x) = (x,7t(x)) for some x £ L^. 

Remark 2.1 For any u £ R and x G U with 7r(x) = v the preimage satisfies 7r~^(v) = Lx- 

The unperturbed equation: We are interested in the ergodic behavior of a strong solution 

of a Levy driven SDE with discontinuous components which takes values in M and respects the 
foliation. In order to avoid that jump increments lead to an exit from the foliation of the initial 
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condition, the jump increments must respect the curved structure of the local coordinates and 
therefore necessarily satisfy a canonical Marcus equation, which are equivalent to the generalized 
Stratonovich equation in the sense of Kurtz, Pardoux and Protter [H]. We consider the formal 
canonical Marcus stochastic differential equation 

dXt = Fo{Xt)dt + F{Xt)odZt + G{Xt)odBt, Xq = xq G M, (4) 

which consists of the following components. 


1. Let Z = {Zt)t^o with Zt = {Z^,... , Zl) be a Levy process with values in for fixed r G N 
on a given hltered probability space Ct = with characteristic triplet (0,z^, 0). 

Suppose that the hltration satisfies the “usual” conditions in the sense of Protter [50] , 

As a consequence of the Levy-Ito decomposition Z is a pure jump process with respect to a 
fj-finite measure v : .8(R”) —>■ [0,oo] called the Levy measures satisfying 



A ll^lP) i'{dz) < oo 


and i^({0}) = 0. 


For details we refer to the monographs of Sato [52] or Applebaum |2]. 


(5) 


2. Let Fq G C'^{M,TTl) with Fq{x) G T^L^- The vector field F G L{W]TTl)) satisfies 

that the map M 3 x ^ T’(x) is and the linear map F{x) sends a vector 2: G R” (-)• F{x)z G 
TxLx to the tangent space of the respective leaf. 


3. 


Let B = {B^,... ,B^) be a standard Brownian motion with values in R” dehned on f2 and 
G G C2(M, L(R”, TIOT)) with G{x) G for any x G M. 


We further assume that the vector fields Fq, F, {DFq)Fq, {DF)F, G and {DG)G are globally 
Lipschitz continuous with Lipschitz constant £ > 0. 

A strong solution of the formal equation (|4|) is defined as a map X : [0, 00 ) x 14 —)• M satisfying 
P-almost surely for all t ^ 0 

w = xo + r Fo{X,)ds + f G{Xs)dBs + ^ [\DG{Xs))G{Xs)d{B)s 
Jo Jo ^ Jo 

+ f FiX,_)dZ,+ ^ ( 6 ) 

0<s^t 

where (B). stands for the quadratic variation process of B in R” and the function <h^^(x) = 
Y (1, x; Fz) and Y (t, x; Fz) for the solution of the ordinary differential equation 

^Y{a) = F(Y{a))z, y(0) = x G M, z G RL (7) 


The perturbed equation: This article studies the situation where an SDE in the sense of ([6|) 
which is invariant on the leaf of the initial condition xq is perturbed by a transversal smooth vector 
field eKdt and the stochastic differentials eG o dB and eK o dZt with e > 0 in the limit of e 0. 
More precisely, we denote by X^, e > 0 the solution in the sense of equation Q of the formal 
perturbed system 

dXf = Fo{Xt)dt + F{Xt) o dZt + G{Xt) o dBt + e(^K{Xf)dt + K{tt{X^)) o dZt + G(7r(Xf)) o dBt ), 

Xq = xo, ( 8 ) 

where the additional coefficients are defined as follows. 
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4. Let iL : M —)• TM be a smooth vector field. 

5. Let Z = (Z^,..., Z’’) be a pure jump Levy process with values in M'’ defined on fl and Levy 
measure i^' satisfying 



z\\‘^)u'{dz) < oo 


and i^^({0}) = 0 


and K eC^iV,L{m\TM)). 

6. Let B = {B ^,..., B^) be an M^-valued Brownian motion defined on and G G C^(y, LiW,TM)). 

We assume that the vector fields K, {DK)K, K, {DK)K, G and {DG)G are globally Lipschitz 
continuous with Lipschitz constant £ > 0. 


Theorem 2.1 (im, Theorem 3.2 and 5.1) Under the preceding setup in particular items a), 
b) and 1.- 3., there is a unique semimartingale X whieh is a strong global solution of ^ on ft in 
the sense of equation It has a eddldg version and is a (strong) Markov proeess. 

Remark 2.2 It is obvious that under the preeeding setup in particular items a), b) and 1.-6., there 
is also a unique strong solution on ft of equation (0) in the analogous sense of equation 
and with the same properties, if Fq is replaced by Fq + eiL and F by {F,£K), G by {G,eG), B by 
{B,B) and Z by {Z,Z) aceordingly. 

We state the crucial chain rule for the Marcus equation given in m (Proposition 4.2). 

Proposition 2.3 Let Z and Fq,F satisfy items 1) and 2) and X be the solution of (0) with initial 
eondition xq with G = 0. Then for any ^ G C^(M'^) we have P-a.s. for all t ^ 0 

{D^){Xs)Fo{Xs)ds + f {D^){Xs-)F{Xs-)odZs. 

Jo 

A direct consequence of the chain rule is the following support property given as Proposition 4.3 
in m- Each jump increment of the noise AZ is mapped to an increment of the solution X 

of ([H). The increment follows the integral curve T in ([7|) along the vector field Fz which is 

tangent to the (smooth) manifold Lxq and hence M. A standard support theorem for ODEs applied 
in the proof of this result then yields that the solution after the jump once again is an element of 
Lxq and hence M. This support property is maintained for an additional Stratonovich component, 
the reasoning is standard. For the solution X‘^,e > 0 given in ([8|) this remains obviously true only 
for M. Under the aforementioned conditions these lines of thought lead to the following foliated 
structure of X: xq G M implies Xt(xo) G Lxg, P-a.s. for all t ^ 0. We shall call a solution of an 
SDE of the type (jH) which admits a foliated solution a foliated Levy diffusion. 

In addition, we obtain that xq G M and e > 0 imply Xf{xo) G M, P-a.s. for all t ^ 0. 


4'(w)=4'(xo) + r 
Jo 


2.2 The main result 


Hypothesis 1: Compactness and Integrability. (a) Any leaf Lxq G 9IT, xq G M, is compact 
and the map xq i —>• diamLxg is Lipschitz continuous in the embedding space of M. 

(b) There is a constant p ^ 2 such that the Levy measures u (of Z) and F (of Z) satisfy 



z\\^ iy{dz) < oo 


and 



z\\'^^ u'{dz) < oo. 
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Hypothesis 2: Existence of invariant measures on each leaf. (a) The solution X of (111) 

has for any initial condition xq G M a unique invariant measure fixo with supp(/ia;(,) = Lxq- 

(b) For vq = 7r(3:o) being the vertical coordinate of some xq G M we define for h : M —)• T5Ji 
with h{x) G TxM 

Q’^{vq) := f h{u)iJ,xoidu) ( 9 ) 

and suppose that for any such function which is globally Lipschitz continuous the 
function D 1/ 9 u e->• Q^(v) G is globally Lipschitz continuous. 

Remark 2.4 Note that (just as L^q) only depends on the vertical component Tr{xo). 

Hypothesis 2 ensures that for each xq G M, vq = 7r(xo) G V the stochastic differential equation 

dw = {w) dt + K{w) o dZt + G{w) o dBt, rc(0) = vq gV (10) 

has a unique strong solution w = i’o))te[o,Too) '^oo being the first exit time of w from V. 

Hypothesis 3: Ergodicity in terms of L^. Let Hypotheses 1 and 2 be satisfied for some p ^ 2. 
We assume that there exists a bounded, continuous, decreasing function p : [0, oo) —)• [0, oo) 
with r](t) \ 0 as t oo such that for any xq G M 

^ r]{t), for all t ^ 0. (11) 


E 


y / TrK{Xs{xo))ds - Q^^{7r{xo)) 

t Jo 


On the rate of convergence. Results about rates of convergence go back to Pascal [48] and Kol¬ 
mogorov [39|, see also [5] and references therein. Recent developments for Levy driven dissipative 
systems can be found in Kulik [40|, see also lailllllBlIlOl. In [401120] for instance the authors 
develop generic methods to establish exponential convergence to an ergodic limit measure in terms 
of the total variation distance for the solution of an SDE driven by a pure jump Markov process. In 
this context the exponential rate of convergence in total variation implies that the rate of conver¬ 
gence in m is of order l/t^. In Subsection 12.31 we provide the simple example of a Levy process on 
the unit circle for which we calculate the precise rate of convergence, which is of the same type. For 
Brownian diffusion processes Hormander’s hypoellipticity condition ensures exponential rates of 
convergence in total variation, see for instance [S] [H [27] and references therein. On the other hand, 
there is no standard rate of convergence for general Markovian systems in the ergodic theorem, see 
for instance Krengel [38] or Kakutani and Petersen [30]. Therefore, it is natural to formulate the 
result in terms of the function 77 following the approach in Freidlin and Wentzell |21j . 

The main result of the article is proved in Section 5 and reads as follows. 

Theorem 2.2 Let Hypotheses 1, 2 and 3 being satisfied for some p 2. Then we have for any 
A G (0,1) and xq G M positive constants Eq £ (0,1), C > 0 and c > 0 such that for any e G (0, Eq] 
and T G [0,1] 


/ 


\ P 


E 

sup \7r{Xl {xo))-w{t)\P 

^CT 

+ ri (cT lne|) 

V 

tG[0,T] ^ 

) 



where is the solution ((3) and w the solution of /f7^) . and = 5^ A Too. 5^ is the first exit time 
of X^{xo) from U in a) and Too is the first exit time of w from V. 


6 








Remark 2.5 Since our main result focuses on Levy processes with only p-th moments, the coeffi¬ 
cient G can he set to 0 since no additional difficulty to the proof of f25\/ shows up. The coefficient G 
will be also dropped in the proof. Including it in the proofs of the sectionslM cmdlSiis straight-forward. 

Remark 2.6 In the proofs of Section 3, 4 cLnd 5, it will turn out that under the preceding assump¬ 
tions none of the constants depends on the precise shape ofV and hence U. Hence without loss of 
generality and for the sake of readability we may assume in the proofs that y = 


2.3 Example: Perturbed Levy rotations of the unit circle 

We illustrate this phenomenon in M = \ {0} with the 1-dimension horizontal circular foliation 

of M where the leaf passing through a point xq € M is given by the (nondegenerate) horizontal 
circle 

Lxo = {(Ikoll cosO, ||xo|| sin6»), 9 e [0, 27r)}. 

Let the process Z = be any pure jump Levy process with second moments. The Levy-Ito 

decomposition of Z yields almost surely for any t ^ 0 


= 



0 ./LKi 


zN{dsdz) -|- 



0 J\z\>l 


zN{dsdz) 


(13) 


where N is the random Poisson measure with intensity measure dt0 n and N denotes its compen¬ 
sated counterpart. Consider the foliated linear SDE on M consisting of random rotations: 


dXt = kXt o dZt, Xq = xq, with A 



(14) 


Equation (|14jl is defined as follows. Note that for some jump increment of Z, z £ M, z 0, we 
have to consider the solution flow 4* of the equation 

-^Y (a) = F{Y {a))z, Y (0) = (x, y), where F{x, y) = A(x, yf', 
obtained by a simple calculation as 


4>^^(x,y) = Y{l-,{x,y)) 


xcos(z) — ysin( 2 ;) \ 
X sin(z)-|-y cos( 2 :) y 


such that 


Xt — xq-\- 


r AXs-zN{dsdz) + Y, - Xs- 

0<s^t 


F{Xs-)AsZ). 


The chain rule of the Marcus integral, Proposition 4.2 in m. states for ||(x, y)^|P := x^ -|- 


d\\Xtf = -2Xt-AXt- o dZt = 0. 


(15) 


In fact, X can be equally defined as the projection of Z on the unit circle. If we identify the plane 
where X takes its values with the complex plane C we obtain Xt = By the Levy-Chinchine 
representation of the characteristic function of Z we obtain for any y G R 

= exp(f'I'(p)), where 'I'(p) = f {e^^^ — 1 — izpl{\z\ ^ l})i>{dz). 
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The invariant measures Hxq in the leaves passing through points xq G M are therefore given by 
normalized Lebesgue measures in the circle centered in 0 with radius ||xo||. We are interested 
in the effective behavior of a small transversal perturbation of order e: 


dXl = AXf o dZt + eK{Xf) dt + edZt 


with initial condition xq 
satisfying 


(1,0), where Z is a pure jump Levy process with Levy measure 



\z\\^u'{dz) < oo. 


We shall consider two classes of perturbing vector fields K. 

(A) Constant perturbation sK = s{Ki,K 2 ) G This example was carried out in [28] for the 
Gamma process on the unit sphere, with Z = 0. The case of a general Levy process is virtually 
identical. The main result in this case reads as follows. For any Levy process Z with E[|Zi|^’] < oo, 
p ^ 2 and A G (0,1) we obtain eq G (0,1) and Tq > 0 such that for any T G [0,To] and e G (0,eo] 
we have 

^ e^T. 


E [ sup 

TTriXf) - 1 

1 

\se[o,T] 

€ 

A 


(B) General linear perturbation eK{x,y) = £A{x,y)'^ = e{ax + by,cx + for a given matrix 
A G which is obviously globally Lipschitz continuous and smooth. The radial component of 

the vector field K is then given by 

'KrK{9, r) = r{{asm{6) + 6cos(0), csin(0) + dcos{6))'^, (sin(0), cos(0))'^) 

= r(asin^(0) + dcos‘^{9) + (5 + c) sin(0) cos(0)), 


where 9 is the angular coordinate of (x, y) whose distance to the origin is r. Hence the average of 
this component with respect to the invariant uniform measure on the leaves (circles) is given by 


1 n -i- d 

r) = —J^ UrK{9, r)d9 = -^r 

for leaves with radius r. We verify the convergence (Hip of Hypothesis 2 for the radial component 
and p = 2. Let Z be a Levy process in with finite fourth moment. Elementary but lengthy 
calculations which can be found in Appendix 16.11 show that 


E 


I/" TrrK{Xs)ds - Q^^^{xo) 
t Jo 


2 >-00 


0 , 


(16) 


where the rate of convergence rj is of order I/'/i as t oo. 

For an initial value xq = (ro cos(0o), sin(0o)) the transversal system stated in Theorem 12.21 
is then w{t) = roe“^*ro. Hence the result guarantees that the radial part 

accelerated time scale ^ has a local behavior close to the exponential in the sense that for any 
A G (0,1) there are constants C,c\ > 0 and Sq G (0,1) such that for any T G [0,1] and e G (0, eo] 
we have 


( E 

sup 

7r^(X|^.^e(xo)) 

2 

V 

sG[0,T] 

e 



^ Cr(e^ + (c| lne|)-5), 


where cx is given in Corollary 13.21 This averaging error tends to zero for fixed T when e 0 and 
for fixed e if T 0. 















3 The perturbation error 

In order to prove the main theorem we have to control the error — X in terms of L^. This result 
relies on the following elementary but in this context crucial lemma on dynamical systems, which 
yields on the right-hand side only quadratic dependence on the “jump increment” z. Due to its 
importance for this article we provide a sketch of proof. 

Lemma 3.1 For a globally Lipschitz continuous matrix-valued vector field F £ L(]R”, 

and z G denote by {Y(t‘, x, Fz))t^o the unique global strong solution of the ordinary differential 
equation 

dY 

— = F{Y)z y(0, X, Fz)=xe (17) 

1) Then there is a eonstant C > 0 such that for any z G M” and x,y £ M with Y{t]x) = 
Y{t;x,Fz) we have 

sup\{DF{Y{t;x))z)F{Y{t;x))z — {DF{Y{t;y))z)F{Y{t;y))z\ ^ C \x — y\ || 2 ;||^. 

2) For any X £ M we have sup^gp,!] \\DF{Y{t; x))F{Y{t] x))\\ < oo. 

Proof: We lighten notation and omit the parameter Fz in Y and write Fz = F^. By the change 
of variables we have for any x £ M, t 0 that 

FfiY{t;x)) = Ffix) + [ DF,{Y{s-x))F,{Y{s-x))ds. 

Jo 

Differentiating in t yields 

j^FfiY{t-x)) = DFfiY{t-x))F,{Y{t-x)). 

Hence for any x,y £ L the mean value theorem and equation (I17p yield 

DFfiY{f,x))FfiY{f,x)) - DFfiY{t-y))F,{Y{t-y)) 

Y d 

= J -^DF^{Y{t]x + a{y - x))){y - x)da 

{D'^Ffi){Y (t; x + a{y - x))){{y - x),F^{Y (t; x + aiy - x)))da. 

Since {DF)F is Lipschitz continuous and F £ the operator D{{DF)F) is uniformly bounded. 
The chain rule D{{DF)F) = {D^F)F + {DF){DF) yields 

\\{D^F)F\\^ ^ \\iD^F)F + (DF)(DF)||oo + ||(DF)(DF)||oo < oo, 

where ||DT|| is uniformly bounded since F is globally Lipschitz continuous and F £ C^. Therefore 

\DFfiY{t;x))F,{Y{t;x)) - DFfiY{t-y))F,iY{t;y))\ ^ \\{D^F)F\\\\zf\x-y\^C\x-y\\\z\\\ 

Since the right-hand side is independent of t we take the supremum as claimed in statement 1). 
Statement 2) is a straight-forward consequence of the product rule. ■ 
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Proposition 3.1 Let the assumptions of Subsection 2.1 and Hypotheses 1, 2 and 3 he satisfied for 
some p 2. Then for any Lipschitz function /i : M —)• M there exist positive constants eo,ko,ki, k 2 
with ko < 1 such that for allT ^ 0 satisfying sqT ^ ko, e £ (0,eo] implies 



sup |/i(Xf (xo)) - h{Xt{xQ))f’ 
.t^T 


^ kiE^ exp(A:2r). 


(18) 


In addition, the constant k 2 is a polynomial in diamLxQ of order p with positive coefficients. 


Corollary 3.2 Let the assumptions of Proposition {Til be satisfied for some p 2. Then for any 
A G (0,1) given there exist positive constants c\, Sq, k^ such that T^ := —cx ln(e), s G (0, eo] satisfies 



sup \h{Xf{xo)) - h{Xt{xo))\P 


< k^e^. 


(19) 


In addition, the constant k^ is a polynomial in diamLxQ of order p with positive coefficients. 


Proof: Plugging T^ = —cln(e) into the right-hand side of (1181) we obtain kisex.p{k 2 T^) = kis^ . 
Given A G (0,1) we choose c = -^ (l — A') and A' = ^ (A -|- l) to infer the desired result. ■ 


Proof: (of Proposition 13.ip The proof consists in three parts. After changing the coordinates 
in part 1 we estimate the transversal component — u| using Lemma 13.11 in part 2. In part 
3 we estimate the horizontal component \u^ — tt| before concluding with a nonlinear comparison 
principle. Part 2 and part 3 are given in separate lemmas. The main tools to derive two (nonlinear) 
comparison principles are Lemma l3. II and Kunita’s maximal inequality for the L^ norm [p ^ 2) of 
the supremum of compensated Poisson integrals found in [42] and an extension of this result for 
p G [1,2] by Saint Loubert Bie [55] . 


I. Change of coordinates: First we rewrite the respective solutions of equation ([1]) and ([8]), 

X and X^, in terms of the coordinates given by the diffeomorphism ip 

{ut,vt) := ip{Xt) and {ufivD := ip{Xf), e G (0,1), t G [0,T]. 

The Lipschitz regularities of h and (p yield a joint Lipschitz constant Cq := Lip{h o ip~^) such that 

IH^t) = \hoip-^{ufivf) - hoip-^{ut,vt)\ 

^ Co\{uf - ut,vf - vt)\ ^ Co{\uf - ut\ + \vf - vt\). (20) 

The proof of the statement consists in calculating estimates for each summand on the right hand 
side of equation above. We define 

S' := (Dip) o F o ip~^, ^ := {Dip) o K o ip~^, A := {Dip) o K o ip~^, 

whose derivatives are uniformly bounded. Considering the components in the image of ip we have: 

^={^H,^v)i ^={^H,^v) 

with G TLxq with Rh T and £ TV ~ with T The chain rule for 

canonical Marcus equations (Theorem 4.2 of |41j ) yields for equation ([8]) the following form of the 
components in ip coordinates 

dul = ^Q{ul,vl)dt + ^{ufivf) o dZt + e^H{ul,vl)dt + ekH{vl) o dZt with uf G L^q, (21) 
dvl = s^v{ufivf)dt + ekv{vf) o dZt with uf G V. (22) 

Note that for e = 0 the equation yields vt = v^ = 0 G V almost surely. However we will write vt 
nevertheless for the sake of readability. 
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Lemma 3.3 (II. Estimate of the transversal deviation \v’: — v. |) Under the previous assump¬ 
tions we obtain the following. There is a constant > 0 such that for eqT < 1 e G (0, Sq] implies 

E[ sup \vf - vt\P] ^ CieP{l + T^P+i). (23) 

t€[ 0 ,T] 

In addition Ci = Ci{diamLxg) depends globally Lipschitz continuously on diamU^Q. 

Proof: By assumption F, K and K are globally Lipschitz continuous. Without loss of generality, 
we can assume they have all a common Lipschitz constant The compactness of and the 
embedding in yields the existence of Ci = diamLxg such that |ttf(xo) — Us{xo)\ < Ci. We use 
the notation 


C 2 {xo) = sup ||jlv(j/,0)|| ^ P(xo)|| + sup l\xQ -y\^ P(a;o)|| +idiamLxo < oo, (24) 

which is finite by the compactness of Lxq and the continuity of ^y. Keeping in mind that 
{Dg{x),u) = p\x\P~‘^{x,u) for x i—>■ g{x) := \x\p,x G we apply the change of variable for¬ 

mula and obtain 


\vt-vt\P=pl \vl-Vs\P -Vs,eKviul,v^,))\ds 
Jo 


+ p f -Vs-,ekv{vl_)odZs)\ 

Jo 

\vl - Vs\P~^\e^v{u%v^^) - £^v{us,Vs)\ds 

Jo 

+ P {vl - Vs\P~^\eSivius,Vs)\ds 

Jo 

+ p [ - Vs-\P~^\{vy - Vs-,£yv{vt-) - kvivs-))dZs)\ 

Jo 

+ p [ - Vs-\P~‘^\{v';_ - Vs-,£kvivs-)dZs)\ 

Jo 

+ p \vl_- Vs-\P-^\^^'^^^^^{vy) - 

0<s^t 


- {vl_ - Vs-) - £{kv{vl_) - kviVs-))AsZ\ 

+ pY I<- - - Vs- - £kv{Vs-)AsZ\ 

0<s^t 

= Hi + H2 + H3 + H^ + H5 + He. 


(Hi) 

m 

m 

{Ha) 

{H,) 

{Ho) 

(25) 


1. Pathwise representation and estimates: Hi: The compactness of yields for 6*3 = 

pl{l + diamLxo) 

Hi ^ £p£ f ^ eCs f \vl - Vs\Pds + eCs f \Vs-Vsf~^ds. 

Jo Jo Jo 

(26) 

H 2 : A direct computation gives 


H 2 ^ 


£pC2 f — Vs\P ^ds ^ eCi f juf — Vs\P ^ds. 

Jo Jo 


(27) 
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H3: Switching to the Poisson random measure representation with respect to the compensated N', 
for instance in Kunita [l2], we obtain for C 5 = pi /|u||>i ll'2||^^(c^-2) 


H3=p[ [ 

Jo Jk 

+ p 

>0 ^||z||>l 

^ ep 






vi_-vs-r\vi 


^{Vs_ — Vs-,s{^v{Vs-) — Av{vs-))z)N'(dsdz) 
— VsY’~'^{vl — Vs,e{k.v{vl) — hv{vs))z)v' {dz)ds 

- Vs-, {kv{vl_) - 'kv{vs-))z)\N'{dsdz) 

+ 6(^5 [ \vl-Vs\^ds. 

Jo 

H4 : With the help of Holder’s inequality we obtain for Cq = pi 

Hi=P f f \Vs- — Vs-\^~‘^{Vs- — Vs-,£^v{vs-)z)N'{dsdz) 

Jo JK’’ 

'fl. 

Jo J U >1 


(28) 


+ P 


Vg — Vs\^ ^{Vg — Vg,eKv{vg)z)u'{dz)ds 


^ ep 



I- -Vg-\^ - Vg-,kv{vg-)z)\N'{dsdz) 


+ eCe / \vl — Vg\^ ^ds. 


(29) 


H5 : For the canonical Marcus terms, Lemma EH statement 1), provides a positive constant C7 
which depends on the leaf of the initial condition such that independent of 0 G [0,1] 

\{Dkv{^^'^"{y,9))z)kv{^^'^"{y,O))z]-\{Dkv{^^^"{x,e))z)kv{^^^"{x,0))z] ^ C 7 \x-y\\\zf. 

(30) 

The Poisson random measure representation of the random sum reads as the following estimate in 
terms of the quadratic variation of Z for Cg = C^pj^ 


0 <s^t 


- Vg- 


m-1 




-{Dkv{^^'^'^^^^(yg))^gZ)ky{^^~^'^^^^{vg))^gZ 
^ e^Cs ^ \vl_-Vg-\^AgZf. 


0<s^t 


The representation of this sum in terms of the Poisson random measure, for instance in Kunita 
], is given for Cg = i)u||>i \\z\\^v'{dz) by 


\vi_-vg-nj^gZf = 


0<5^t 

which yields 



\vl_ — Vg-J’\\z\\^N'{dsdz) + Cij f \Vg—Vg\^ds, 

Jo 


H5 ^ £2^8 


// 

Jo JK’ 


vl_ — Vg-\^\\z\fN{dsdz) + e"^Cio / \vl — Vg\^ds. 


/‘ 


(31) 


12 










He : For the last term we use Lemma l3.11 statement 2), which yields a positive constant Ci\ such 
that for any z G M*' 


sup \\{Dhv{Y{0,{),£hvz))zAv{Y{9,^,£hvz))z\\ < Cii|| 2 :|p < oo. 

0 G[O,1] 

Hence exploiting that J||^||>i < oo we have for C 12 = C'ii(p/2) /|| 2 ||>i 

He S^P - Vs- - £'kv{Vs-)^sZ\ 


0<s^t 

2 

^ ^ sup ||(My(OA,Z)^v(OA,Z)|| 
z 

0 G [ O ,1] 


^=Y{0,Vs,eSiv^BZ) Qc^s^t 




= b‘Cu 


0 <s^t 

fl 

Jo JM’ 


V,_ - Vs- 


1^ ^\\z\\^N'{dsdz) + £‘^Ci 3 f 

Jo 


— UoF ^ds. 


Combining the estimates ([26l [271 [28l [29l [STJ [32]) we obtain 


\vt - Vtl"^ ^ £€3 f \vl — Vs\^ds + £63 [ + f \vl - Vs\^ ^ds 

Jo Jo Jo 


(32) 


(33) 


+ £Ce [ \vs - Vs\^ds + £Ce [ 

Jo Jo 


— Us|^ ^ds 


+ e^Cio [ |uf - Us|^ ds + f - Vgf ^ds 

10 Jo 

|uf_ - - Vs-, (kv{vl_) - kv{vs-))z)\N'{dsdz) (34) 


+ £p 


fl 

Jo Jm.^ 



+ £p I I \Vs--Vs-f ^\{vl--Vs-,kv{vs-)z)\N'{dsdz) 
rt 


+ e^Cs 


+ £^<712 


\vl_-Vs-\P\\zfN'{dsdz) 


fl 

Jo Jr" 


\vl_ - Vs-\P-f\zfN'{dsdz). 


(35) 

(36) 

(37) 


2. Estimate of the marginal expectation: Taking the expectation we obtain a constant C 14 
snch that 


iE[i< - vtin ^ ^^14 r (iE[i< - ^.r] +iE[i 

Jo 


Vs - Vs\ 


ip-l 


])ds 


and 


_ p-i. 


sup E[|uf - Us|P] ^ eCu / (E[|uf - + E[|uf - p )ds 

sG[0,t] ^0 

^ £Ciit sup E[|uf - Usl^] + Cii£ f E[|uf - ds. 

se[o, 4 ] Jo 


, P-1 
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For e G (0,eo] and t G [0,T], T > 0 such that £oTCi 4 ^ that is fixing /cq = ( 2 C 14 ) ^ in the 
statement, we have 

sup lE[|uf — ^ 2(7146 f sup E[|uf — vt\^]~ds. 

se[o,r] Jo te[o,s] 

It is easy to verify that the maximal solution of this equation is given by 

sup E[|uf - ^ Ci 5 (er)P. (38) 

sg[o,t] ^ P ' 

We can replace the exponent p in the estimate of E[|uf — by 2p and p — 1. This is possible 
since the integrals with respect to v' do not depend on p and obtain for q G {2p,p — 1} 

sup E[|u| - ^ (7i6(er)'^. (39) 

s&[0,T] 


3. Estimate of the expectation of the supremum: We go back to (I33p and note that all 

integrands with respect to the Lebesgue integral ds are positive, such that the integrals are positive 
and increasing. Further we note that the last four summands are compensated Poisson random 
integrals M = + M^. The sum of the terms ()341 [35l [Mil will be denoted for convenience by 

and term (|37l) will be denoted by We obtain with the help of Jensen’s inequality 

E[sup — u|^] ^ eCi 7 / (E[|uf — — Vs\^]~)ds + E[sup 2 + E[sup 

[o,T] Jo [o,T] [o,r] 

(40) 


For we use Kunita’s maximal inequality for the exponent 2 (see [l2] or m) and Young’s 
inequality for the exponents and 2p in order to obtain 

E[sup |M^|^] ^ 6(7i8(||.^y(0)|| / [ E[\vl - Vs\'^^^~^'^]\\z\\'^i''{dz)ds 

[o,r] ^ Jo JK'- 

+ r [ E[\vl-vsnzfu'{dz)ds 
Jo Jk’" 


+ / E[\v^^ - Vs\^^]\\z\\^iy {dz)ds 

Jo JR’’ 

T 


^ eCi 9 f 

Jo 

^ C'lo J ^(1 + ^)IE 




p-i) 


+ E 


V, - V, 


|2p 


ds 


Ts -Vsl 


|2p 


+ e^^'jds, 


and inserting (l39|) in (iTT|) we get 

E[sup |m1|2] ^ C 20 [^{esf^ds + ^ (72oe^^r(l + T^^- 

\ 0 .T] Jo 


(41) 


(42) 


This could be also repeated with term but the price to pay would be f || 2 ;||®z 2 '((i 2 ;) < 00 , 

equivalent to the finiteness eighth moments of Z. Instead for we use a maximal inequality 

for integrals with respect to Poisson random measures for the exponent 1 given in |49] . Lemma 
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8.22 and resp. Theorem 8.23, going back to Saint Loubert Bie [53], which states the existence of a 
constant C 21 such that for C 22 = <^21 and t ^ 0 


E[sup|m2|] 

[o,r] 




£^<^22 


fl 

Jo Jw 


E 


V, - V, 




\z\\^i 2 '{dz)ds < £^6*22 


j-T 


/ E 

\vl-Vsf’ 

Jo 



n 

^ ds. 


Inserting (f38l) in the preceding expression we get 

E[sup |m2|] ^ e^C 23 riesf-^ds ^ 

[0,T] Jo 


(43) 


Transforming (I40p with the help of (I38p . (|42p and (1431) and keeping in mind that eT < 1 yields 

rT 


E[sup \v^ - ^ £^24 / ((es)P + {esy-^)ds + C 23 eP^/T{l + T^ + 

[o,r] Jo 

^ (7266^(1 TT^P+i). (44) 


Lemma 3.4 (III. Estimate of the horizontal component \v!: — u.\) Under the previous as¬ 
sumptions we obtain the following. There are constants C2,to > 0 such that for e^T < 1 e € (0,eo] 
implies 


E 


sup |ttf — Utf 
tG[0,T] 


< C2 £^p{1 + (T V exp (Car). 


(45) 


Proof: For convenience of notation we restart with the numbering of constants. Formally we 
obtain 


u, -Ut = 


Ut= idoiul,Vs)-doius,Vs))ds+ / {d{ul_,v‘;_)-d{us-,Vs-))odZs 

Jo Jo 

+ £ [ - ^Hius,Vs))ds + e j JiH{us,Vs)ds + e j kH{vl_)odZs. (46) 


10 

This equation is defined in M” as 
rt 


ul-ut = 


[ [do{u%v^^)-do{us,Vs)]ds 

Jo 

+ j [d{ul_,v';_) - d{us-,Vs-)]dZs 

+ ^ [(<I.S^=^«_,uf_)-#^^^(u,_,u,_)) 


0 <S^t 


- «_ -Us-,vl_ - Vs-) - (?«_,<_) -d{Us-,Vs-))AsZ] 

+ £ j {^h{u%vI) - ^H{Us,Vs))ds + £ I {Us, Vs)ds + £ kH{vl_)dZs 
Jo Jo Jo 

0<s^t 
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The change of variable formula for (I46p yields formally 


\ul-ut\P=p / \ul-Us\P - Us,do{ul,vl) - ^o{Us,Vs))ds 

Jo 

+ p [ \ul_ - Us-\P~‘^{ul_ - Us-,{d{ul_,vl_) - d{Us-,Vs-)) O dZs) 

Jo 

+ ep \ul-Us\^~‘^{ul-Us,AH{ul,vl)-JiH{us,Vs))ds 

Jo 

+ ep \ul - Us\^~‘^{ul - Us,^HiUs,Vs))ds 

Jo 

+ P \ul_ - Us-\^~‘^{ul_ - Us-,ekHivl_) O dZs). 


This is defined in M” as 

rt 


ul-ut\^=p \ul-Us\^ ‘^{ul - Us,So{ul,vt) - ^o{'<J-s,Vs))ds 

Jo 

+ p [ \ul_ - Us-\^~‘^{ul_ - Us-,idiul_,vl_) - dius-,Vs-))dZs) 

Jo 

0<s^t 

- «_ - Us-,vl_ - Vs-) - -d{Us-,Vs-))AsZ) 

+ ep \Us - Us\^~'^{Us - Us,JiH{Us,vl) - JiH{Us,Vs))ds 

Jo 

+ sp \ul - Us\^~'^{Ug - Us,'^HiUs,Vs))ds 


+ ep [ \ul_-Us-f ^{ul_ - Us-,kHivl_)dZs) 

Jo 


+ pY 1 

0 <s^t 

- ivl_ - Vs-) - siknivs-) - kHivs-))AsZ) 

+ p ^ \ul_ - Us-f~‘^{ul_ - Us-,^''^^^‘’^iVs-) - Vs- - ekHivs-)AsZ) 


ih) 

ih) 

ih) 

ih) 

ih) 

ih) 

{h) 

ih) 


0 <s^t 

=■ Ii +12 + I 3 + h + I 5 + h + Ij + h- (47) 

We now estimate the eight summands on the right-hand side one by one. The estimates Ii and 
I 4 are straight forward Lipschitz estimates. For the stochastic ltd terms we use the different kinds 
of maximal inequalities. The estimate of the canonical Marcus terms I 3 , /y and /§ is the most 
laborious task in which we exploit the result of Lemma 13.11 The term is straight forward. 

1. Estimate of the stochastic Ito integral terms I 2 and Jg: I 2 : Due to the existence of 

moments of order at least 1, I 2 has the following representation with respect to the compensated 
Poisson random measure associated to Z 

Ul_ - Us-Y’~‘^{ul_ - Us-, - ^{Us-,Vs-))dZs) 
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+ 


{ul_ - Us-, (^(u^_,v^_) - S(us-,Vs-))z)N(dsdz) 

n K - Usl^~^(ul - Us, (S(Us,v^) -^(us,Vs))z)u(dz)ds. 
z\\>l 


(48) 

(49) 


For the first term (j48p we apply the embedding L? <Z L^, Kunita’s maximal inequality (see [l2] or 
m) for exponent equal to 2, the compactness of L^q implying |M^(a:o) — n. (xo)| ^ diamLxQ, and the 
elementary Young inequality for the exponents p and pjip — 1) combined with inequality (1381) . We 
obtain a positive constant Ci > 0 such that for = 2C'i(/jgr WzW^i'{dz)y^'^diarriLxo and C 3 = C 23 
we have 


E 


sup 

te[o,r] 


^ CiE 
^ CiE 


/ / \ul_-Us-\^ '^{ul_-Us-,{d{ul-,Vs-)-d{us-,Vs-))z)N{dsdz) 

Jo JR’’ 

/ / \ul-Us\‘^^P~‘^^\{ul-Us,{^{ul,vl)-d{us,Vs))z)\‘^v{dz)ds 

Jo Jr^ 

[ [ \ul - Us\‘^^^~^U\ul - Usl"^ + \vl - Vs\‘^')\\z\\‘^u{dz)ds 

Jo 4r’' ^ ^ 

1 

dsV +C3ePTP+^. 




( / 

sup \u^ — u\P 

V Jo 

L [0,.] J 


(50) 


The second term is less delicate. Young’s inequality for the exponents p/{p — 1) and p and yields 


E 


sup 

t£[0,T]J0 4||^||>1 


- Us\P ^(u^ - Us, {^{ul,vl) - ^{us,Vs))z)v{dz)ds 




||z||z/((iz)E sup / (\Ug — Us\P + \Ug — u, 
|| 2 ||>i '-te[o,r]4o 




^ £ WzWv 
4|bll>i 

^ (74 / E 

Jo [ 0 ,s] 


/ r'^ 


r'^ 


2 / E 

sup \u^ — u\P 

ds+ E 

\VI-Vs\^ 

^ Jo 

L [O.s] 

Jo 



\)d. 

ds 


sup — u\P 


ds + C 5 ePTP+\ 


(51) 


le: Let N' be the compensated Poisson random measure associated to the Levy process Z. Then 


sup e / \ul_-Us-\P ‘^{ul_ - Us-,^H{vl-)dZs) 
£fo,ri Jo 


iG[o,r] 

= sup e f f \ul_ - Us-\P~'^{ul_ - Us-,{kH{vl_) - kH{vs-))z)N'{dsdz) 
tG[0,T] Jo Jr'" 

[ [ \ul-Us\P~'^{u^s-U s,{^h{v^s) - ^H{vs))z)u'{dz)ds 

Jo J\\z 


+ sup e 

te[0,T] Jo 4||^||>1 


Ug_ — Us-\P ^{Ug_ — Us-,^h{vs-)z)N'{ dsdz) 


+ sup e 

tG[ 0 ,r] Jo 4 R’’ 

+ sup e / \ul — Us\P~'^{Ug — Us,hH{vs)z)i''{dz)ds. 
tG[o,r] Jo 7|b||>i 


(^ 1 ) 

iJ2) 

(Jo) 
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The terms Ji and J 2 are structurally identical to (jdHjl and and are estimated analogously to 
(150]) and (I5TI1 where 5 is replaced by which yield the estimates 


E 


sup \£ 
te[o,T] Jo 


f! 

Jo Jw 


\ul_ - Us-\^ - Us-, {kniVs-) - kH{vs-))z)N'{dsdz)\ 




E 


10 


sup — u\^ 


[o,d 


ds\ + 2 and 


E 



sup e| / / \ul-Us\^ Us,{kH{vl)-kH{vs))z)i''{dz)ds\ 

'-te[o,T] JO J|h||>i 
T 


/ E 


/ 


sup — u\^ 


[o,d 


ds + C9£PTP+\ 


For the term J 3 we observe that = 0 such that ^v{vs) is constant. Kunita’s maximal inequality 
for the exponent 2 yields the constant Cio and the boundedness of — Us| the constant Cn = 

C'lo (/iR»- \\z\\'^v'{dz){diamLxQ)P~‘^^'^ such that 

j-t 

E 


sup £\ 
t&[ 0 ,T] 


^ eE 


sup 

te[o,T] 


^ sCio 
^ eCii 


[ [ \ul_-Us-\P '^{ul_-Us-,kH{vs-)z)N'{dsdz)\ 

Jo JK’' 

/ [ \ul_ - Us-\P~'^{ul_ - Us-,kH{Vs-)z)N'{dsdz)f 

Jo Jk.’" 

'\ul - Us\‘^^P-^A\\z\\^i^'{dz)ds 


E 


sup \u^ — u\P 

[0,s] 


ds 


The 


term J 4 is again easier, for C 12 = i|| 2 ||>i we obtain 


E[ sup £ 


t&[0,T] Jo J||2||>1 

j-T 

/O [0,s] 

rT 


ul - Us\P ^{ul - Us,^H{vs)z)v'{dz)ds\ 




£Ci 2 [ E[sup 

JO [0,s] 

C 12 / E[sup \u^ - u\P]ds + Ci 2 £PT. 
Jo fo.sl 


This yields 


E[sup I/el] ^ Cis 

[0,T] 


/ E 

/o '-tG[0,s] 


sup |u^ — u\P 


*) 


E 


tG[0,s] 


sup \u^-u\P ds + £PTP{2 + T) 


(52) 


2. Estimate of the canonical Marcns terms I-j and /g: Lemma 13.11 tells us that 

there is a positive constant, ( 714 , say, which depends on the leaf of the initial condition such that 
independent of 0 G [ 0 , 1 ] 

{D:S{^^\y,e))z):S{^^\y,e))z] - 0))^! ^ Cm\x - y\\\z\f. (53) 
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Now we apply the Taylor’s theorem to 9 i—)• Y{6]x,^A.sZ) = {x,9) with 0 G [0,1] and 

inequality (1531) . The Poisson random measure representation of the random sum writes as the 
following estimate in terms of the quadratic variation of Z 

0<s^t 

- K_ -Us-,vl_ - Vs-) - idiul_,vl_) -5K-,^^s-))As2')| 


^2 ^ {Dm^^^^K-,yU))^sZ)d{^^^^^{ui_,vU)^sZ 

0<s^t 

- {Dd{^^^^^{us -, n,_)) AsZ)^{^^^^^{us -, Vs-)) AsZ 

^C'i4( ^ \ul_ - Us-l^WAsZW^ + ^ \ul_ - Us-\^~^\v^s_ - Vs-\\\AsZ\\ 
0<s^t 0<s^t 

^2Cu( Y1 \<-- ^s-HAsZf + Y, \vl_-Vs-nAsZ\ 


0<s^i 


0<s^t 


(54) 


The representation of this sum in terms of the Poisson random measure, for instance in Kunita 
I, of the first term is 


Y Wl-- Us-Y’\\AsZ\ 


0<s^t 


/7 

Jo Jw 


u^_ — Us-r z 


P||.l|2 


N{dsdz) + f f 

Jo J\\: 


ul — ns|^||z||^ v{dz) ds 


(55) 


0 J\\z\\>l 


and the analogous result if |n^_ — u^-j is replaced by |nf_ — n<j_|. The maximal inequality for 
integrals with respect to the compensated Poisson random measures and inequality (|38p yield 


E[sup j/sj] ^ Cis f f ^EPsup — nj^j + E[|nf — Us]^]') ||z||^ z^((iz) ds 
[o,T] Jo Jr^ ^ [0,s] 2 

= <^15 [ lk|pi^(d-2)( [ (e[sup — ttj^] + E[|?;^ — ds 

Jr'" ^ Jo ^ [o,sl 2 


(56) 


^C'lef f E[sup jrt^ — rtjP] ds + . 

^ Jo [0,s] 2 

I 7 : For If we use Lemma l3. II statement 1) in terms of ()30p and Young’s inequality 
^ |n^_ - Us-\P-H<- - Us-,<^^'^^^^^{vU) - ^^'^«^^^{vs-) 

0<s^t 

- {Vs- - Vs-) - sikniVs-) - kH{Vs-))AsZ) 

^ Y \<- - ns-r>^^^^^HvU) - - «_ - Vs-) - e(kH{vl-) - kH{vs-))AsZ\ 




0<S^£ 
2. 




\ul_-Us-\^ Vf- - 


0<S^t 


0<s^t 
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We rewrite the last expression in terms of the compensated Poisson random measure N' and obtain 


E 


0<s^t 


III J' 

i:l 


■u,_|p)||A,Z ||2 


^l_-Vs-\P)\\zfN'{dsdz) 

(57) 

+ — Vs\P) \\z\\‘^h''{dz)ds. 

(58) 


+ 

/o ^lhl|>i 

The maximal inequality in [53] by Saint Loubert Bie for the exponent 1 yields 


E 


/7 

Jo Jw 


sup 
te[o,T] Jo 


“ W-T + 7 s- “ 'f^s-D \\z\\^N'{dsdz) 


^ Cis 


fl 

Jo Jm.^ 


E 


- Usf + \vl - 


z\\‘^i''{dz)ds 


^ Cig f WzW^V^dz) 

Jw 

■f 


E 


sup \u^ — u\^ 

[0,s] 


ds+ E 


'4-Vs\P 


ds 


^ C20 / E 
/o 


sup \u^ — u\^ 

4G[0,s] 


ds + C 2 iePTP+\ 


The term (j58]l is treated obviously such that 

E[sup I/ 7 I] ^ £^ 6*22 [ 

[o,T] Jk’' Jo 


sup — u\P 


[0,s] 


ds + C23£^TP+\ 


Ig: For Ig Lemma [3Tl statement 2), yields 

^ |<_ - Us-\P~^{u^s- - - Vs- - ekH{Vs-)^sZ) 

0<sCi 

^ K--Us-?~^\^^'^^^^Hvs-)-Vs--eAH{Vs-)/J^sZ\ 

^£2^24 ^ \ul_ - Us-\P-^\\^sZ\\\ 


0 <s^t 


(59) 


such that Saint Loubert Bie’s maximal inequality with exponent 1 and elementary Young’s estimate 
for parameters p yield 


[0,T] 



rT 

5 / W^fj^'idz) / 

jk’’ 

Jo 



/ E 

sup \u^ — u\P 

Jo 

L [o,s] J 


sup \u^ — u\P ^ 

[0,s] 

ds + C 2 q£PT . 


ds 


( 60 ) 


3. Estimate of 

i-T 


f 


\ul-Us\P - Us,£^H{Us,Vs))ds ^ C27 / e\ul-Us\P ^ds 

Jo 

^C27 [ £\ul - Us\Pds + C27£PT 

Jo 
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such that 


E[sup I/ 5 I] ^ £(727 [ E 

[0,T] Jo 


sup \u'^ — u\^ 


L [0,s] 


ds + C 27 SPT. 


(61) 


4. Nonlinear comparison principle: Taking the supremum and the expectation of the left- 
hand side of equation ()171) and combining the estimates of ®"[sup[o,T] 1-^*1] given by (l50l) . (f^ . 
(|5^ . (l59|) . (| 6 fl and (fHTIi we obtain a positive constant C 28 


E 

- 1 

e 

1 

^C'28f 

f\ 

sup \u'^ — u\P 

ds ( 


sup \u^ — u\P 


te[o,r] 

V 

Jo 

te[o,s] 

V 

Jo 

t€[ 0 ,s] 


For the concave invertible function G{x) = x + \/x, := E 

ePfP^l 1 ) we have achieved the integral inequality 

rt 


supte[o,r] Wt - ut\^ 


ds ] + ePTP{l + T) ). 

and p{t) = 


V’(i) ^ C2sG(^ j 'ijj{s)ds'j + G28P{t). 

Hence the convex function G~^{x) = (x -|- ^) — x — \ yields for C 29 = max{2C28,4(7|g} 
G-\m) ^ G-^{c 28 G[j\{s)ds) +G 28 P{t) 


^ ]^G-^{ 2 G 28 G(^j\{s)ds)^+]^G-\ 2 G 28 P{t)) 

^G29 f ^is)ds + G29G-^{p{t)). 

Jo 


Note that G ^{x) ^ kx for all n G (0,1) and x ^ xq for some xq > 0, which we calculate as follows. 


1 . 


1 


{xo + o) - \/ ^0 + 7 = KXO Xo = 


(1 — kY 


The quadratic function qx‘^ satisfying that qx‘f^ = kxq has the prefactor 

qx^ = kxq q = — = (1 — k)^, 

Xq 

and we obtain that any k G (0,1) and x > 0 satisfy G~^{x) ^ min{(l — tdf'x, kx}. Hence 

min{(l - ^ (7“^(T/'(t)) ^ (729 f 'ilj{s)ds + C29G~^ {p{t)). 

Jo 

2 

This implies for 'ip{t) ^ kxq = the integral inequality 

'ip{t) ^ — [ ilj{s)ds + — 

^ Jto ^ 

for t ^ to, where to = inf{t > 0 | V'(t) = kxq}. Gronwall’s inequality yields for t ^ to 

, G29' 


'ilj{t)^^G ^(p(t)) exp (^(t - to)^). 


21 
















Since the running supremum 11 —>■ ^(t) is monotonically increasing we also obtain that for any t ^ 0 

-ipit) < — {p{t) V p{to)) exp t] . 

Ki \ tv / 

Taking into account that G~^{x) ^ for all x ^ 0 we obtain for ^ and Cso = 2 C 29 


E 


sup \u^ — u^’ 
[0,T] 


^ C' 3 oe'^(l + (T V exp (C 30 T). 


(62) 


IV. End of the proof of Proposition 13.11 Eventually Minkowski’s inequality, the Lipschitz 
estimate (l20l) . the sum of the vertical (1231) and the horizontal (HSl) estimate yield the desired result 
for ^3 = Cl V C 2 and constants ki,k 2 > ^ 


E 


sup \h{Xl{xo)) - h{Xt{xQ))\P 
iG[0,T] 


^ C3 + T^P+i) + C3oe"^(l + (T V exp (C3or)) 

^ /cie^exp {k2T). 


This finishes the proof. ■ 

4 The averaging error 

For convenience we fix the following notation. Given /i : M —)• M” a globally Lipschitz continuous 
function and : V —)> M” its average on the leaves defined as ([9]). For f ^ 0, xq G M and e G (0,1) 
denote the error term 

■= [ /i(V|(xo)) - Q^(7r(X|(xo)))ds. 

Jo ® ® 

Proposition 4.1 Let the assumptions of Provosition \3.1\ he satisfied for a fixed p 2. Then for 
any h : M ^ ML globally Lipschitz continuous, A G (0,1) and xq G M there are constants hi > 0 
and So G (0,1) such that for e G (0,eo] T G [0,1] we have 


/ 




E 

sup \6^^{e,s)\P 

) 4biT 

+ p (cT lne|) 

V 

sE[0,T] 

) 



where c = c\^p is given by Corollary \3.2\ and p is the ergodic rate of convergence given in equation 
071) of Hypothesis 3. 

Proof of Proposition (|4.1|) : The common part : Fix xq G M. For e G (0,1) and T > 0 we 

define the partition 

T 

to=Q<t\<-- - <t%e 

with the following step size 

Ag := -cT ln(e), 
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where c > 0 is given by Corollary 13.21 The grid points are defined as 


:= nAg for 0 ^ n ^ N^, e € ( 0 , 1 ) and the total number is 
We rewrite the hrst summand of 6 ^ by 


Ne=[ 


1 


ce\ ln(e)| 


J+1. 


f h{Xf{xo))ds = e [ h{X^{xo))ds 
Jo ^ Jo 


^ ftn+l 

= e^ / h{Xl{xo))ds + e h{Xl{xQ))ds. 

We lighten notation and omit for convenience in the sequel the superscript e and h as well as the 
initial value xq whenever possible. The triangle inequality yields 

^ + I^2(T^)I + I^3(T^)I) (63) 


where 


A^£ — 1 , 

Ai{t,e) j [/i(Af(xo)) - h(X,_t„(Xf^(xo)))] ds, 

n =0 

A^e —1 ,1 

/ ‘' 71+1 

A 2 {t,e) := e X] / IH^s-tAXlixo))) - Q( 7 r(Xf^( xq)))] ds, 

n =0 

d^ 3 ii,£) ■■= y] £^eQ{'^{Xt^ixo))) - / < 3 ( 7 r(X|(xo))) ds. 

n =0 do 

The following three lemmas estimate the preceding terms. This being done the proof of Proposi¬ 
tion [Q is hnished. 


Lemma 4.2 For any A G (0,1) and xq G M there are 62 > 0 and Eq G (0,1) such that for any 
e G (0, eo] and T G [0,1] 


( E 

sup | + (s,e)|^’ 

V 

se[o,r] 


1 

P 


^ b2T£^. 


Proof: Fix A G (0,1). We apply the Markov property of Xt and Holder’s inequality in the time 
variable and the fact that by dehnition ^ T 




p Ng — l 


P^n+1 1 - 

E 

sup |+(s,e)|P 

= e 

E 

1 / [h(Xf(xo)) - h{Xs-tAXlM))] ds\^ 1 


se[o,r] 

n =0 


J tn 


r r 

= e^EE\ [hiXl_,Jy)) - h{Xs-tAy))] ds\^ I y = Afjxo) 

n =0 ^ ^ do 


< eNs A^E 


E 


sup \h{XI{Xl{xo))) - h{X,{Xl{xo)))f 

sG[0,ti] 


y = ^LM 


1 

p 
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^ T max E 

nG{l,...,Afe} 


E 


sup \h{Xl{Xl^{xQ))) - h{Xs{Xl^{xQ)))\^ 

sG[0,ti] 


^ T max 
nG{l,...,Afe} 


¥\kz{diam{LxUxo)))]£^ ■ 


y = ^Li^o) 


The last estimate is an application of Corollarv l3.2l where we have found that the value = k 3 {dyg), 
dyg = diam^Lyg) from (fT^ is a polynomial in of order p with positive coefficients. 

First we estimate the term involving k^. Since by Hypothesis 1 the mapping vq i—>■ dy^ is Lipschitz 
continuous there is a Lipschitz constant I' on 17, which yields the constants ai and 02 such that 
almost surely 


hidx^J ^ ks{dya + £' dist{LxQ, Lxejxo))) 

^ ai{k3{dyg) + 1 ) + a2dist{Lxo,Lx^jxo)T- ( 64 ) 


The next step consists in the estimate of the term E[dist(L 2 ;p, (xq))^]- Using Xt{xo) G Lxq 
for all t ^ 0 Lemma EJ ensures the existence of the positive constant C'i(d„Q) which is an affine 
function in dy^ such that for k = C'oC'i(d,;o) 


E[ sup dist{Lxo, Lxe(^xo)f]^ ^ ^qE 

t 6 [ 0 ,ti] 


sup |uf(xo) - Us(xo)|^ 
se[o,ti] 


< ke^ 


(65) 


where Cq is the Lipschitz constant of the local coordinates (p with Lipschitz constant 7". We obtain 

E[ sup dist(L„o,Lx^e(,^j,))P]F ^ E[ sup dist(L^o, + E[E[ sup dist{Ly^, Lxi(xo)T \ ■ 

tG[o,t2] te[o,ti] te[ti,t2] 


The hrst term on the right-hand side obeys (|65l) . the second one can be calculated recursively 
E[E[ sup dist{Lyo,Lxe^xo)T I ^ E[E[ sup dist{Lxo, Lxe^y)T \ y = X^^{xo)]]p 

t€[ti,t2] iG[0,4i] 

^ CoE[Ci{dxi_^{xo))]^^ 

^ {k + fke^)e^. 


Hence 


E[ sup 

tS [ 0 ,t 2 ] 


dist{Ly^, Lxi{xo)T]^ ^ 


We argue by induction 


E[ sup dist{Ly^,Lx^(xo)T]^ 
telo,t„+i] 

^E[ sup dist{Lyg,Lx^[xo)f]^ -hE[E[ sup dist {Ly^, Lx^(y)f \ y = XfJxo)]]p 

< M( ^ (() GoE[<;i{rfv..(»))]£". 

^ i=0 ^ ^ 

We continue with the second term 

CoE[Ci(djfe (3,Q))]e^ ^ C'oCi(d„o)e^-F7"E[ sup dist{Ly^, Lx^i^xo))]^^ 

t^[0,tn] 
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+ sup dist{L^Q,Lx^(xo)f]^^^ 

tS [0,tTl] 


< —i 
^ £" 


i=0 




eventually leading to 


,1 k 


E[_sup d\si{Lx,,Lxiix,)Y]~^ ^ -^)+y>[Yl 


tG[0,tjVe + l] 


Ne 


i=0 
Ne + 1 


N, 


i" 


Ne + l 


i=l 


Ne 

i + 1 


{feY" - 1 




for £ G (0, eo] sufficiently small. Combining the preceding result with (IM|) yields for e G (0, eo] 


E 


sup |^i(s,e)|P 
se[o,r] 


^ T max W\k^{diam{Lxe 

^T£^(ai{k 3 {dyo)+ 1 ) + 02^1 sup dist{Lxo, Lxi(xo)T]) 


^ T£^(^ai{k3{dyf^) + 1 ) + a2{kN^£^y^, 

and a positive constant 62 which yields the desired result for any e G (0, eo] and A G (0, A — p(l — A)) 


E 


sup |^i(s,e)|P 
se[o,T] 


< b2T£^. 


Lemma 4.3 Let t] be the rate of convergence defined in m- For the process A 2 in inequality LU^) . 
T > 0 fixed and A G (0,1) here is c = c(A) G (0,1) for which there are Eq G (0,1) and 63 > 0 such 
that for any T G [0,1] and e G (0, eo] we have 


( E 

sup 1^2(5, e)r 

V 

se[o,T] 


1 

P 


< 63T g (cT| lne|). 


The proof is virtually identical to the proof of Lemma 3.3 in [25j in the purely Brownian case. 


Lemma 4.4 Due to Hypothesis 2 is globally Lipschitz continuous. Then the process A 3 in 
inequality satisfies the following. For any A G (0,1) there are eo G (0,1) and 64 > 0 such that 
for any e G (0, eo] and T G [0,1] 


( E 

sup 1^3(5, e)|^’ 

V 

se[o,r] 


4 

P 


^ bfiTe^. 
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Proof: We calculate 


1 ^' 


^£-1 rNeAe 

(r,e)| = eA,Q(7r(XfJ) - / g(7r(X|)) ds 

' ZZ, Jo 


n=0 
Ne-1 

sup \Q{tt{XI)) - Q{7r{Xl))\ 

^_Q £tn^S<C.£tn-\-l 

Ne-1 

^eA,cY^ sup |uf-ufj. 


( 66 ) 


Minkowski’s inequality, the Markov property, Lemma 13.31 with Ci in d^Q and estimate (1651) lead to 
E[ sup |.43(s,e)|^’]r ^ TC max E[E[ sup IvUtAv) - Voiy)\^ \ V = Nni^oW 

/V-J- + . 


5G[0,T] 


nG{l,...,Ne} tn^s<tn+i 

^TC max E[E[sup \v';{y) - v^o{y)\P \ y = Xt„{xo)]]p 
ne{l,...,iV£} 


^ TCE 


Ci{diam{Lx,^ (xo))) 


-Ns 

^TCCi{d,,){e^ + m,£^^) 

^ Tbi£^ 

for A G (0,2A-1). ■ 

(Proof of Proposition 14.11 Combining Minkowski’s inequality with Lemma 14.2114.41 yields for 
any A G (0,1) and p ^ 2 constants CA,p, eoi^o S (0,1) and bi > 0 such that for any T G [0,1] 
satisfying eqT < ko e £ (0, Sq] implies 


E 


sup |<5^o(e,i)r 

t&[0,T] 


^ 6ir(e^ + r/(cr|ln(e)|). 


5 Proof of the main result 

We keep the notation of the Proof of Proposition 14.11 By the change of variable formula for 
canonical Marcus integrals, m Proposition 4.2, we may rewrite 


v^t = 


[ Kv{u\_,u%_)ds + f {kv){v%_) o dZg. 

^ Jo e e Jq e 

Since equation (fT0]l tells us that w[t) = (w(s))ds + Jq hv{w{s—)) odZg we obtain 

v\ - w(t) 

e 

= [ {w{s))ds + f {Av{v%_) — Av{w{s—))) o dZs 

io ® ® Jo ^ 

= ^xir(£)^)+ [ Q^^iu\,vi) - Q^^{w{s))ds + f [kvivl_) - kv{w{s-)))dZs 
Jo ^ ^ Jo ^ 

+ _ (^^(^1^) _ Kviw{s-)))Aj^ 

0<s^t 
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— ^xo (E) i) + Ol +02 + 03- 


Now, for T E [0,1] we obtain with the help of Jensen’s inequality 


E 


sup|Oi|P ^aiTP-^ 
[0,T] J 


E 


sup \v\ — w{t)\P 

iG[0,s] ^ 


ds. 


Kunita’s maximal inequality [32] yields a positive constant 02 = 02 (p) such that for the constant 
^3 = 02 ^(/]Rr \\z\\‘^i''{dzy/^ + 2 \\z\\Piy'{dz)) we have 


E 


sup \O 2 \P 
[0,T] 


+ E 


sup 
te[o,T] Jo 


+ E 


< a 2 E 


sup 

tG[0,T] Jo ./|h||>l 
■T 


+ a 2 lE 


(^\/(u|_) — Av{w{s—)))zN'{ dsdz)\P 

[kv{v%) — hv {w{s))) zu' {dz)ds\P 

[ f l^kvivl) — kviw{s)))z\Pi''{dz)ds 
Jo Jm.^ ® 

ill \{kvivi) - kviw{s)))z\'^iy'{dz)ds] 

''{dz)E f \kv{v%) — kv{w{s))\^ds 
L Jo ^ 


lhll>i 




/ E 

sup |ui — rc(s) ^ 

Jo 

■sG[0,t] ® 


dt. 


Finally the Lipschitz continuity of the vector fields ky and {Dky)ky and Lemma l+T] provide a 
constant 04 > 0 


lOsr ^ I Y. 


0<s^T 




- (u|- - w{s-)) - {eky{v^^) - £ky{w{s-)))AsZ 


I V \iDkyivU)AsZ)kyivU)AsZ - {Dkyiwis-))AsZ)kyiwis-))AsZ\y 

Z _ e e / 


o<ss£r 


^(04 Y - ^(s-)iii^s^ii' 

o<s^r 

^ (a4)P( Y “ ^(s“)l 11^.5^11 

0<s^T 


(67) 


Switching to the representation in terms of the Poisson random measure 


Y. -'P^{s-)\\AsZy 

0<s^T 




e 


w{s—)\\\zy N\dsdz) 



tc(s)|||z|p v'{dz) ds, 


( 68 ) 
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we obtain 


]E[sup| 03 r] ^ 2P-^(a4f E 

[ 0 ,T] V 


sup 

te[o,T] 



vl- — w{s—)\\\z\\‘^N'{dsdz)\^ 


+ E 




e 


n;(s)|||zf 



—: O 4 + O 5 . 


( 69 ) 


We apply Kunita’s maximal inequality [2] which yields a constant 05 = 05 ( 1 ?) > 0 and Jensen’s 
inequality. For = 2 ^’“^(a 4 a 5 )P and 07 = aeifj^r + (fj^r hence obtain 






sup |n|—t(;(s)|^ 

sG[0,t] ® 


dt. 


(70) 


The term O 5 follows straight forward. Summing up X)i=i E[sup^g[Q 7 ^] \Oi\^] we obtain a constant 
ag such that for T e [0,1] 




rT 


sup \v\ 

-w{t)r 

^ ag / E 

sup |ni — t(;(s)|P 

^t&lO,T] ® 


Jo 

_6G[0,t] " 


dt + E 




The standard nonautonomous version of Gronwall’s lemma implies a constant 09 > 0 such that for 
T G [0,1] 


E 


sup |u|_ — r(;(s—)|^ 
se[o,r] ® 


^ E 


6^nT,Er] expiagT) ^ awE\d^^{T,sr 


Finally an application of Proposition 14.II for h = Ky finishes the proof of Theorem 12.21 


6 Appendix 

6.1 Detailed calculations of the example 

In the sequel we verify (fT 6 l) . Keeping in mind Q^''^{xo) = we obtain 


E 


1 /■* 

- / 'KrK{Xs)di 

^ Jo 

r [ 

- / E (asin^(Zs) + (icos^(Zs) + (6 + c) sin(Zs) cos(Zs)) 
t Jo 

pi 

-at-\ — / E ((d — a) cos^(Zs) + (6 + c) sin(Zs) cos(Zs)) 
t t Jq I 


ds 

ds 


( 71 ) 
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= ra + 


jid-a) f 

r Jo 

/' 


E 


cos^{Zs) ds + {b + c) / E[sin(Zs) cos(Zs)]ds 
-I Jo 


= ra + -{d — a) 




a + d d — a r 

- r ^ - 

2 2 wo 

a + d d — a r 
- r H- 

2 2 i ./o 


exp(i 2 Zs) 

r /"* 

- / 3?exp(—s'I’(2)) (Js 
^ Jo 


d — o 

as H-- —r 


f exp(—s3ft'I'(2)) cos(s3ft'I'(2)) ds. 
Jo 


Since 3?'I'(2) > 0 we have 


E 


1 

- / 'KrK{Zs)d. 

^ Jo 


a + d 
s - r 




|d — a| r 


^ J exp(—s'I'(2)) ds 


Id — o| r 
2^'(2) t 


0, t ^ CK 


This shows the result for p = 1. For p = 2 we continue 


E 

= E 


1 /■* 

- / TrrK{Zs)d. 

* Jo 


a + d 2 i 
s - r 


= r^E 


= E 


f r(asin^(Zs) + dcos^{Zs) + (6 + c) sin(Zs) cos(Zs))ds- — — r 

Jo 2 J 

J (asin^(Zs) + dcos^(Zs) + (5 + c) sin(Zs) cos(Zs))ds^ 

— (a + d)- f (asin^(Zs) + dcos^(Zs) + (6 + c) sin(Zs) cos(Z 5 ))ds + ("———') 

t Jo V 2 / 

1 /■* \ 2 - 

- / (asin^(Zs) + dcos^(Zs) + (6 + c) sin(Zs) cos(Zs))ds j 

LVt Jo / J 

(a + d) E - J (asin^(Zs) + dcos^(Zs) + (6 + c) sin(Zs) cos(Zs))ds ^ 



t^OO by (dH) 


We calculate directly 


Et =E 


J (asin^(Zs) + dcos^(Zs) + (5 + c) sin(Zs) cos(Zs))ds^ 


2 i 




t 2 



a^E 


'0 ao 


sin^(Zs) sin2(Zo-) 


;„2/ 


+ d^E 


COS^(Zs) COS^(Zo-) 


+ (6 + c)^E sin(Zs) cos(Zs) sin(Zo-) cos(Zo-) 


+ adE 


sin^(Zs) cos^(Zo-)] + adE cos^(Zs) sin^(Zo-) 


+ a(6 + c)E sin^(Zs) sin(Zo-) cos(Zo-) + a(6 + c)E sin(Zs) cos(Zs) sin^(Zo-) 
+ d(6 + c)E cos^(Zs) sin(Zo-) cos(Zo-) + d(5 + c)E sin(Zs) cos(Zs) cos^(Zo-) 
We apply the elementary identities 

cos(x)^ cos^(y) = ^ ( cos(2(x — y)) + cos(2(x + y)) + 2 cos(2x) + 2 cos(2y) + 2^ 


dads. 
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sin(x)^ sin^(y) = ^ cos(2(x — y)) + cos(2(x + ?/)) — 2 cos(23:) — 2 cos(2y) + 2^ 

sin(x)^ cos^(y) = ( cos(2(x — y)) + cos(2(x + y)) + 2 cos(2x) — 2 cos(2y) — 2^ 

sin(x)^ sin(y) cos(y) = g ( sm(2(x — y)) + sin(2(x + y)) + 2 sin(?/) + 2^ 

cos(x)^ sm{y) cos(y) = “ g ( sin(2(a; — y)) + sin(2(x + y)) — 2 sm{y) + 2^ 

and obtain 


1 /■* f r 

Jo Jo 8 


Et = 




(f ( 

+ -{« 

, jb + cf 

-{[m 


^i2{Zs—Z(j) 

+ m 

^i2[Zs-\-Zo-) 

-25RE 

^i2Zs 

- 25RE[e*2^" 

'^i2{Zs-Za)' 

+ 5RE 

^i2[Zs^Za) 

+ 2$RE 

g*2Z/ 

+ 2$RE 

g2Za 


5RE 


^i2{Zs-Z,y) 


-m 


^i2{Zs-^Z(j) 


J2{Zs-Z„) 


+ 5RE 


^i2(Zs+Za) 


+ 2m 


^2Zs 


1 a(6 + c) 

9E 

gi2(Z«-Z^)" 

-9=E 

^i2{Zs-\-Za-) 

+ 29E 

4 

d{b + c) 

4 1 

L 

9E 

gi2(Zs-Z<,)' 

- 9E 

^i2{Zs-\-Za) 

- 29E 


- 25RE 

g2Z, I 

J2Za^ ' 


J2Za 


+ 2 


- 2 


Since 


E 


gi2(Zs±Zo-) 

= E 

gi2Zs±CT 

^ g-(6±a)xl>(2) 

and 

E 

gi4Zs 


dads. 


= g-s'^(4) 


for s ^ cj we get 
rt rt „2 


Et = 


1 r r r 

“2 / / — j exp(—|s — (7|5R'I'(2)) cos(|s — (t| 9'I'(2)) + exp(—(s + cj)$R'I'(2)) cos((s + (j)$R'I'(2)) 

E Jo Jo 8 1 


+ 


— 2exp(—s5?^(2)) cos(s9‘'I'(2)) — 2exp(— (T$R'I'(2)) cos((79'I'(2)) + 2| 

— \ exp(—15 — (j|5R^^(2)) cos(|s — (j|S>^(2)) + exp (—(5 + a)3?^(2)) cos((5 + a)3^(2)) 

8 I 

+ 2exp(—s5R'I'(2)) cos(sSj'I'(2)) + 2exp(—( t5R'I'(2)) cos(it9''I'(2)) + 2| 

— exp(—|s — c7|5R'I'(2)) cos(|s — C7|9''I'(2)) + exp(—(s + cj)'b(2)) cos((s + (t)9‘'I'(2)) 

+ 2exp(—s5?'I'(2)) cos(s9''I'(2)) — 2exp(—(T$R'I'(2)) cos(s9''I'(2)) — 2| 

{b + cf 
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a{b + c) 


+ 


+ 


+ 


I exp(—|s — (t|$R'I'(2)) cos(|s — (t| 9'I'(2)) — exp(—(s + (7)3ft^(2)) cos((s + cj) 9'I'(2))| 
I exp(—|s — (t| 5R'I'(2)) sin(— |s — (t| 9''I'(2)) + exp((s + cr)3?'I'(2)) sin((s + cr)9''I'(2)) 


— 2exp(—fT'I'(2)) sin((T9‘ 
d{b + c 


— 2 exp(— 


I exp(—|s — cj|$R'b(2)) sin(|s — cr|9'I'(2)) — exp(—(s + cr)3ft'I'(2)) sin((s + cj) 9'I'(2)) 
p(—cj5R'b(2)) sin(cr$j^'(2))| dads. 
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Since all integrals over the exponential terms converge, their contribution in the preceding sum 
vanishes as t —)• oo as these are divided by and only the constants under the integrals survive. 
Therefore 


such that 
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t—>-oo Cl 
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ad (a + 
Y “ 4 


' 7rrK(Z,)ds - 
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